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FOREWORD 


This  report  was  prepared  under  Expenditure  Order  903“ 
li+2g  by  Captain  Jean  G.  Goppert,  as  project  engineer  in  the 
Wind  Tunnel  Branch,  at  the  request  of  the  Chief  of  the  Project 
Unit,  Wind  Tunnel  Branch. 

Included  among  those  who  cooperated  in  reviewing  and 
editing  the  report  were  Dr.  B.  H.  Goethert  and  Mr.  R.  R.  Kassner 
of  the  Wind  Tunnel  Branch. 
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ABSTRACT 


A  method  is  presented  for  computing  the  lateral  extent  of  a  weak 
shock  wave  due  to  a  small  wedge-shaped  irregularity  of  a  wall  in  two- 
dimensional  supersonic  flow.  Computations  have  been  carried  out  over 
the  Mach  number  range  1.05  to  1.25  for  angular  irregularities  of  the 
wall  of  0.1  to  0,k  degrees.  The  entropy  change  in  the  shock  wave  as 
well  as  the  influence  of  the  wall  boundary  layer  were  assumed  to  be 
negligible.  Computations  were  carried  out  for  the  straight  part  of 
the  shock  wave  and  for  such  points  along  the  curved  part  of  the  shock 
wave  at  which  the  local  Mach  number  change,  A  M,  across  the  shock 
wave,  is  reduced  to  values  of  0.001,.  0.002,  0.003,  0.001+,  and  0.005. 


PUBLICATION  REVIEW 

Manuscript  Copy  of  this  report  has  been  reviewed  and  found  satis¬ 
factory  for  publication. 

FOR  THE  COMMANDING  GENERALs 
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SYMBOLS 


M  Mach  number 

Free  stream  Mach  number 

Mg  (Up  A  Mg)  Mach  number  behind  the  shock  -wave,  but  in  front  of  the 

expansion  fan 

H  (M-i-AMz)  Mach  number  along  some  expansion  wave  at  which  the  flow 
'  has  attained  a  Mach  number  arbitrarily  close  to  Mx. 

£  Length  of  the  wedge 

8  ,/edge  angle 

Ah  Maximum  wedge  thickness 

Q  Angle  measured  counterclockwise  from  the  horizontal 

0,  Mach  angle  of  flow  at  % 

0g  Angle  of  first  expansion  wave 

0-  Angle  of  the  expansion  wave  along  which  the  flow  has  attained 
the  Mach  number 

€  Angle  of  the  shock  wave  at  ary  point 

€«  Angle  of  the  shock  wave  before  intersection  with  the  expansion 

fan 

CL  Mach  angle  measured  from  the  local  flow  direction 

Deflection  angle  required  to  change  the  flow  Mach  number  from 
V\  m  =  1.0  to  M-l  in  a  Meyer,  expansion 

V  Deflection  angle  required  to  change  the  flow  Mach  number  from 
M  =  1*0  to  Mg 

r  Radius  from  the  rear  corner  of  the  wedge 

r2  Radius  to  the  point  of  intersection  of  the  shock  wave  and  the 
first  expansion  wave 

r,  Radius  to  the  point  on  the  shock  wave  where  the  Mach  number  behind 

'  the  shock  has  attained  the  value 
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Lateral  distance  (perpendicular  to  free  stream) 

Distance  from  wall  to  point  of  intersec, si on  of  shock  wave  and 
first  expansion  wave 


Distance  from  wall  to  point  along  shock  wa re  where  the  Mach  number 

behind  the  shock  has  attained  the  value  M, 

? 

C 

_E  Ratio  of  specific  heats 


SECTION  I 


PURPOSE 


The  purpose  of  this  investigation  is  to  formulate  an  estimate  within 
practical  accuracy  of  the  magnitude  and  extent  of  the  disturbance  due  to  a 
small  irregularity  of  a  wall  in  two  dimensional  supersonic  flow.  Considera¬ 
tion  is  limited  to  low  supersonic  Mach  numbers  and  small  irregularities,  ior 
which  only  attached  shock  waves  appear  and  the  shock  losses  are  negligible. 
Boundary  layer  influences  are  not  included,  but  it  can  be  expected  that  t  y 
will  weaken  the  intensity  of  the  flow  disturbance. 


SECTION  n 


DERIVATION 


The  problem  has  been  simplified  by  considering  the  special  case  of  a 
plane  surface  with  a  wedge  shaped  irregularity  in  supersonic  flow. 

Figure  I  shows  the  wedge  shaped  irregularity  AB.  Upstream  of  A  and 
downstream  of  B  the  wall  is  assumed  to  be  plane  and  parallel  to  the  ^ndis“ 
turbed  flow.  A  shock  wave  originates  at  A,  changing  the  flow  angle  from 
to  8  and  the  Mach  number  from  %  to  Mg.  These  flow  conditions,  however, 
exist  only  within  the  triangle  ABPa.  At  point  P2  (wall  distance  Yg)  the 
shock  wave  is  intersected  by  the  first  of  the  expansion  waves  start 

from  the  corner  at  B.  Beyond  P2  the  shock  wave  is  curved  due  to  intersec 
tion  with  the  expansion  waves  and  the  intensity  of  the  shock  wave  ^reduced. 
Therefore,  the  Mach  number  change,  A  M3,  (M-j-  Hj),  across  the  shoe  , 

occurring  at  a  point,  P,  (wall  distance  Y3),  on  the  curved  part,  is  smaller 
than  the  change  A  Mg  across  the  straight  part  of  the  shock. 

The  object  of  the  following  computation  is  to  determine: 


a.  Mach  number  disturbance,  A  M3,  across  the  curved  part  of  the 
shock  as  a  function  of  distance  Y3. 

b.  Mach  number  disturbance  across  the  straight  part  (equal  to 
disturbance  at  wall)  and  distance  Yg. 


Curved  Part  of  the  Shock  wave 

Figure  2  shows  an  element  of  the  shock  wave  between  two  expansion  roves 
having  the  angles  0  and  0  -<j0  with  the  free  stream  direction.  The  shock 
rove  element  forms  the  angle  €  with  the  free  stream  direction. 
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Therefore,  since  (0  -  € )  is  the  angle  between  the  shock  wave  element  and  an 
expansion  wave,  and  r  is  the  radius  from  the  corner,  B,  we  have: 


dr 

-rd0 


cot  (02  “  €  ) 


integrating  from  to  r^  (See  Figure  1) 


(0-€) 


d0 


(1) 


(2) 


Thus  the  first  problem  is  to  set  up  (0  -€  )  in  terms  of  0  .  For  a  given 
initial  Mach  number,  M]_,  and  for  assumed  values  of  the  shock  wave  angle,  €  , 
the  Mach  number,  M,  behind  the  shock  wave  can  be  computed  from  the  equation 
(Eq.  1|.33>  Ref.  3). 


M2 


£1 

1+2 


cos 


YU i 


1  *  ^  b. 


T 


sir?  € 


(3) 


The  proper  range  of  assumed  values  of  €  was  determined  from  a  plot  of  8  vs€ 
for  constant  values  of  taken  from  Table  I  of  Ref.  2.  By  neglecting  the 
entropy  changes  in  the  weak  oblique  shocks,  these  computed  values  of  M, 
together  with  Table  II  of  Ref,  1,  for  issentropic  flow,  can  be  used  to  obtain 
0  ,  which  is  the  angle  of  the  expansion  wave  along  which  the  Mach  number 
is  M.  1’hat  is: 


0  -  CL  *  ( I/j  —  1/ )  .(h) 

where  Ct  is  the  Mach  angle  corresponding  to  the  Mach  number,  M,  and  (]/-!/) 
is  the  flow  inclination  along  the  Mach  wave  with  respect  to  the  direction  of 
the  undisturbed  flow  with  M  =  M^.  (See  figure  3 •) 

Thus  we  have  calculated  values  of  0  and  €  corresponding  to  various 
values  of  M.  It  is  now  possible  to  compute  and  plot  (0  -  € )  vs  0  • 


If  the  plot  of  ( 0  -  € )  vs  0  can  be  approximated  by  a  straight  line  over 
the  desired  range  of  0  ,  then  the  integral  (Eq.  2)  can  be  evaluated  very 
easily  since  we  can  substitute 


9  -  €  ~ 


(5) 


where  the  constants 
as 


a  and  c  can  be  evaluated  from  the  geometry  of  the  line 


a 


. M— 

d  (0  -  €) 
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Similarly,  a  slight  curvature  of  the  plot  ( 5  -  € )  vs  Q  can  be  approxi¬ 
mated  by  two  straight  lines  and  integrated  in  two  intervals.  If,  however, 
the  curvature  is  too  great,  then  a  graphical  integration  seems  to  be  the  best 
procedure.  Merely  plot  a  curve  of  cot  (0  -€)  vs  Q  ,  compute  the  area  under 
the  curve  between  the  limits  Q  2  and  Q-z}  and  take  the  antilog  to  get  £3  . 

r2 

The  angle  Q  2  is  defined  as  the  angle  of  the  first  expansion  wave  of  the 
Prandtl-Meyer  expansion,  (Figures  1  and  3)* 

Thus  Q2  »  a2  +  $  (6) 

Qy.  is  defined  as  the  angle  of  the  expansion  wave  along  which  the  flow  has 
reached  ary  arbitrary  Mach  number,  MV,  which  is  determined  by  the  choice  of 
a  value  of  AM,  across  the  shock  wave.  After  a  selection  of  AM^  is  made, 

6  j  can  be  computed  try  means  of  Eq.  3  and  i+. 


Straight  Part  of  the  Shock  Wave 

r2 

The  next  step  is  to  compute  the  value  of  -j-  which  is  the  radius  from 

the  expansion  corner,  B,  to  the  point  of  the  intersection  of  the  shock  wave 
with  the  first  expansion  wave,  expressed  in  terms  of  the  length  of  the  wedge. 
This  is  a  simple  problem  in  trigonometry  since  the  initial  shock  wave  angle, 

€ 2»  and  the  Mach  number  behind  the  shock,  ML,  can  be  determined  from  tables 
in  Ref.  2,  Eq.  3.  Since  the  direction  of  a  Mach  wave  is  given  by  the  Mach 
angle,  which  depends  upon  the  known  Mach  number  of  the  flow,  we  can,  by 
referring  to  Figure  1,  see  that  the  law  of  sines  will  give  the  desired  radius, 

r2 


T  ! 


sin  (go  - 8  \ 
■sin  (&  2  -€2) 


(2)  e2  -  «2 

(3)  3  =  _ 

X  c 


CL  - 
2  2 


_ sin  (€  2  ***  8  ) 

cos  8  sin  (a2  -€2  +  & ) 


y 

It  is  now  a  simple  step  to  determine  -A  (See  Figure  1) 


h  m  h 
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RESULTS  AND  DISCUSSION 

Computations  were  carried  out  for  the  Mach  number  range  1.05  to  1.25,  and 
for  wedge  angles  of  0.1,  0.2,  O.J,  and  O.lj.  degrees,  using  arbitrary  values  of 

*  .001  to  .005.  fhe  results  are  presented  in  Figures  ^  and  5. 

The  probable  error  in  Yp  and  Yz  may  be  as  great  as  3%  near  Mi  ■  I.05, 
with  considerable  improvement  as  the  Mach  number  increases.  Ihis  large  error 
is  due  to  the  extremely  small  difference  between  the  shock  wave  angle  and  the 
angle  of  the  expansion  waves.  However,  it  must  be  remembered  that  only  the 
approximate  order  of  magnitude  of  the  lateral  extent  of  the  disturbance  is 
all  that  was  to  be  obtained  from  this  calculation. 


The  case  of  negative  values  of  the  wedge  angle,  8  ,  was  checked  at 
^1  “  The  values  of  Y?  were  found  to  be  of  the  same  order  of  magnitude 

but  slightly  smaller  than  tne  values  for  positive  angles.  The  difference  was 
so  small,  however,  that  further  computations  did  not  appear  to  be  warranted. 


The  results  have  been  presented  in  Figure  5  as  a  single  family  of  curves 
of  A  vs  Y  j/x>  8  for  lines  of  constant  Mach  number,  M^»  The  computations 

showed  that  Y^/i  was  proportional  to  8  within  the  investigated  range  of  wedge 
angles  and  Lach  numbers  and  within  the  accuracy  of  these  computations.  It 
should  be  noted,  also,  that  the  values  given  in  Figure  5  are  not  valid  for 

values  of  Y,  <  Y2,  (see  Figure  1),  since  the  disturbance  aM  is  constant 
0  to  Y  ■ i2. 


from  Y 


The  distance,  Yp,  as  well  as  the  disturbance  A  M^,,  valid  for  0  =  Y  =  Yo 
has  been  plotted  in  Figure  h  with  the  wedge  angle  8  as  a  parameter.  For  2 
the  purpose  of  interpolating  between  the  curves,  it  may  be  remarked  that  at 
constant  Mach  number,  %,  AMp  is  approximately  proportional  to  8  ,  while 
Yg/Ah  is  approximately  proportional  to  1/8  •  furthermore,  for  a  given  wedge. 

the  distance  Y2  varies  approximately  with  the  Mach  number  according  to 


APPUCATION  OF  RESULTS 

There  are  two  problems  which  can  be  solved  quite  easily  using  Figure  it. 

a.  Problem  1.  Given  a  wall  configuration,  8  and  1,  and  a  free 
stream  Mach  number,  M-j_.  Find  the  strength  of  the  disturbance,  A  Mz,  at  some 
given  distance,  Y?,  from  the  wall.  3* 

y 

Solution.  Compute  »  and  read  the  value  of  A  from  the 

curve  on  Figure  5  corresponding  to  the  proper  value  of  Mx. 
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b.  Problem  2.  Given  the  maximum  allowable  disturbance,  A Hj ,  at  a 
distance,  Yz,  from  the  wall,  and  a  free  stream  Mach  number,  M, .  Find  the 
allowable  height  of  wall  disturbance  Ah  «  1  8  which  would  satisfy  the  maximum 
allowable  disturbance  requirements. 

Solution,  From  Figure  5>  find  the  value  of  h  corresponding 

to  the  given  values  of  H]_  and  AM^«  ihen,  since  Yj  is  given,  this  determines 
the  allowable  disturbance  height,  X  *  8  • 

If,  in  the  problem  to  be  solved,  the  wall  has  a  curved  irregularity  in¬ 
stead  of  the  assumed  abrupt  wedge-shaped  irregularity,  the  results  given  in 
this  report  will  be  conservative  since  the  initial  shock  wave  will  be  split 
up  into  a  number  of  compression  waves  which  will  probably  converge  into  a 
shock  wave ,  but  the  intensity  will  be  less  than  for  the  wedge. 

The  same  effect  could  be  expected  in  considering  the  boundary  layer, 
since  the  boundary  layer  would  tend  to  round  off  the  corners  of  the  wedge. 


SUMMARY 


The  local  Mach  number  disturbances  in  two  dimensional  flow,  due  to  a 
wedge  shaped  wall  irregularity,  are  calculated  for  various  lateral  distances 
from  the  wall  for  the  Mach  number  range  between  M  *  1,05  and  1.25,  (See 
Figures  1*  and  5«) 

When  the  height  of  the  disturbance  Ah  -  8  *1  is  kept  constant,  the 
lateral  extent  of  the  region  between  the  initial  shock  wave  and  the  first 
expansion  wave  increases  inversely  with  the  square  of  the  wedge  angle  of  the 
wall  irregularity  and  varies  with  the  Mach  number  according  to  _  i  , 

(See  Figure  lu) 

In  the  region  behind  the  first  intersection  of  initial  shock  wave  and 
expansion  wave,  the  local  Mach  number  disturbance  is  nearly  independent  of 
Mach  number  and  wedge  angle ,  but  rather  depends  mainly  on  the  ratio  of  the 
lateral  distances  to  the  disturbance  height,  (See  Figure  5«) 
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Figure  4:  Lateral .distance  to  first  intersection  of  shock  wave  and  expansion  fan  and  Mach  number 
turbance  at  the  wall  for  small  wedge-shaped  irregularity  of  a  two-dimensional  wall. 
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Figure  5:  Lateral  variation  of  Mach  number  disturbance  due  to  a  small  wedge-shaped  irregularity  of  a  two- 
wave  1  Wal1’  Reglon  downstream  of  first  intersection  between  shock  wave  and  expand™ 
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